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SUMMARY 


The purpose of this analysis is to investigate the use of an 
automatic control to damp undesirable lateral oscillations. The type of 
control under consideration is a yaw damper servomechanism with rudder 
control applied in proportion to yawing velocity. The effect of angle 
of attack and gyro axis angle on components of roll and yaw picked up is 
developed, and the optimum angle of tilt of the gyro axis is selected. 

The problems of maintaining equal effectiveness of the control through the 
speed range, eliminating yaw compensation during steady turning or rolling 
maneuvers, critical time lag, and prevention of feedback to the pilot's 
controls are considered. 

The airplane used for the analysis is an operational jet fighter 
of high performance which has been reported to have-a "snaking." or "Dutch 
Roll" problem. The flight conditions investigated are cruising at 35,000 
feet, 504 normal rated power at 10,000 feet, and vower approach at sea 
level. Calculations in these configurations predict that the aircraft has 
an unsatisfactory "Dutch Rolling" tendency, probably in all flight condi- 
tions. Installation of the proposed damper system will enable the airplane 
to satisfactorily meet the lateral-directional oscillation requirements. 

In addition, the stabilizer will cause the aircraft to be spirally stable, 
and will have the effect of gearing the rudder to the ailerons for improved 


entry into rolling or turning maneuvers. 
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SYMBOLS AND COEFFICIENTS 


angle of attack measured from the flight path to the airplane 
reference axis, degrees. 


angle between the reference axis and the horizontal axis, positive 
when the reference axis is above the horizontal axis, degrees 


angle of the flight path to the horizontal, positive in a climb, 
degrees 


angle of attack of the principal longitudinal axis of the airplane, 
positive when the principal axis is above the flight path at the 
nose, degrees 


angle between the reference axis and the principal axis, positive 
when the reference axis is above the principal axis, degrees 


angle between the gyro spin axis and the reference axis, positive 
when the reference axis is above the gyro axis, degrees 


angle between the flight path and the gyro axis, positive when the 
gyro axis is above the flight path, degrees 


angle of aximuth, radians 

angle of bank, radians 

angle of sideslip, radians (v/V) 

sideslip velocity along the Y axis, feet per second 
true airspeed, feet per second 

mass density of air, slugs per cubic feet 

density ratio, p/P, 

dynamic pressure, pounds per square foot (1/2 pV2) 
weight of airplane in pounds 

mass of airplane in slugs, W/g 

—er of gravity, feet per second per second 
relative density factor. m/pSb 

factor for dimensionalization, 2b u/V 


trim lift coefficient, W cos ¥ /q3 


s” 
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mach number 
wing span in feet 


distance from the center of gravity of the airplane to the center 
of pressure of the fin, feet 


wing area in square feet 

control surface area for yaw damper in square feet 
vertical tail area in square feet 

rudder area in square feet 


distance from the longitudinal body axis to the center of pressure 
of the vertical control surface 


the operator, d/dt 
period of oscillation in seconds 


time for amplitude of oscillation to damp to one-half of its 
original amplitude in seconds 


cycles for the amplitude of oscillation to damp to one-half of its 
original amplitude 
° a. r 
control gearing, or gain, | sy | 
ye 
amplitude ratio of damper, 1/K or | /s | Damper 


airplane amplitude, |‘’/ € = 


moment of inertia about the principal longitudinal axis in slug 
feet squared ° 


moment of inertia about the principal normal axis in slug feet squared 


moment of inertia coefficient about the principal longitudinal axis, 


I /qSb 
O 


moment of inertia coefficient about the flight-path axis, 
(I) cos*y + I sin? ) 


moment of inertia coefficient about axis normal to flight path, 
(% cos + I! sin’ ) 


moment of inertia coefficient about the principal normal axis. I /q8b 
O 
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I vroduct of inertia coefficient with respect to flight path axis and 
4 axis normal to flight path, a - I' )sin cos 


g mon 
C, rolling moment coefficient asa) 
C yawing moment coefficient (Yewing moment ) 
, qsb 
Cy side force coefficient ae) 


ry angular yawing velocity in radians per second 


Dy g angular rolling velocity in radians per second 








aCy dc 
Cp > C) = 
B dp B dp 
dCs dc 
Cy = C x n 
re) 7 ag 
aC, dc 
: “np 
P asa) P a(sy) 
aC ac 
Oe _ oa eo 
OLy Cy ay “ny Ca 
= b a b 
he CO Dy “ng * ov 
-_-_ ae 
“Ye ap 1 af 
C 3 a C = a 
Y ag Y pb 
6. r d ev 


The standard NACA conventions in regard to directions of axes, angles, 


control deflections, forces, moments, and linear and angular motions are used. 
2 


Where dotted quantities are written, they indicate time derivatives, ie, = 2, 
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INTRODUCT ION 


The purpose of this report is to conduct a theoretical investigation 
into the use of a yaw damper servomechanism to improve the lateral dynamics 
of an operational high speed jet fighter. Calculations, flight tests, and 
pilots reports indicate that this aircraft does not satisfactorily meet the 
requirements of NAVAER SR-119B, "Specifications For Flying Qualities of 
Piloted Airplanes," in regard to the damping of the lateral-directional os- 
cillation, controls free or fixed. 

High performance requirements of present day aircraft and those of 
the immediate future are leading to designs in which lateral dynamic flight 
stability is becoming a serious problem. "Dutch roll" at low airspeeds, 
particularly on landing approaches, is annoying, fatiguing, and dangerous 
to the pilot. It makes instrument approaches particularly difficult. The 
present day trend toward sweepback and low aspect ratio causes high effective 
dihedral at low apeeds and also reduces the damping in roll and increases the 
adverse yaw due to roll. These factors combine to contribute to unsatis- 
factory lateral dynamics. "Snaking," the control free lateral dynamic oscil- 
lation, compromises the effectiveness of the aircraft in carrying out its 
mission as a gunnery or bombing platform. Satisfactory solutions to the 
problems of "Dutch roll," and "snaking" will not only improve the flying qual- 
ities, but will give the designers more freedom in still higher performance 
designs. 

Considerable work is being done by various investigators on the causes 
of these unsatisfactory lateral oscillations. Aerodynamic balance of the 


rudder is a probable major cause. Reference 2 traced the poor damping qualities 
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of its test airplane (Navy AD) to the rudder hinge moment characteristics. 
Reference 3 points out as a possible cause the effect of fuel sloshing on 
airplane dynamics. Non-linearity of the directional stability and damping- 
in-yaw cerivatives is shown to be a possible cause in reference 4. Random 
separation phenomena such as flow at the tail-fuselage juncture and at the 
tail pipe have been demonstrated to have important effects. Reference 5 
shows that positive inclination of the principal axis can have considerable 
stabilizing influence. Here wing incidence has an important effect. It has 
been demonstrated theoretically that atmospheric turbulence can initiate and 
maintain lateral oscillations. 

The desiger is faced with an imposing problem in selecting the proper 
hinge moment parameters, which will satisfy control force and other require- 
ments, and still not invite snaking. fKeference 2 discusses this problem. 
Irreversible power-boosted aware! systems offer a promising solution. In 
any case, it has become apparent that some means must be found to augment 
aerodynamic damping. The increment needed is considerably in excess of the 
amount that can be provided by an increase of vertical tail area or tail 
length. It can be shown that the tail volume coefficients should be deter- 
mined by the requirements for static sideslip stability, unsymmetrical power, 
and stalling qualities, and that other methods should be sought to add to 
aerodynamic damping. 

Although hinge moment parameters combined with friction and slack in 
the control systems explain conventional rudder snaking in many cases, others 
are on record in which the rudder has been locked or observed not to move. 


Usually the calculated and tested control fixed periods and times to damp 


2 
172 
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agree quite well. The case under consideration in this paver is one of the 
latter group. The manufacturer felt that the lateral oscillations were trace- 
able to non-linearities in the derivatives and to separation effects. It was 
originally thought that fuel sloshing was a contributory cause, but baffling 
was installed and found to have no noticeable effect. In their tests they 
confirmed the effectiveness of a "de-boost" tab, as reported in reference @. 
Since there are no known aerodynamic fixes: for the poor damoing in the rudder 
fixed condition, an automatic stabilization device provides the most promising 
solution. 

It is the purpose of this analysis to consider the effect of using a 
servo stabilizing control, and to select an optimum setting of gain and spin 
axis angle for this application. Two systems are already commercially devel- 
oped (by Minnesapolis-Honeywell and Boeing for the XB-47, reference 8, and by 
Lear for use in conjunction with the F-5 autopilot, reference 6). The theo- 
retical investigation is here conducted with the use of one of these types of 
servomechanisms in mind. 

Reference 7 considers stabilizing devices sensitive to various quan- 
tities and concludes that the only one which has the capability of satisfac- 
torily meeting the requirements of reference 1 is an autopilot which supplies 
rudder control proportional yawing velocity. This is the type of yaw damper 
referred to above. This reference and the above mentioned systems neglect 
the problem of the tilt of the gyro axis. The gyro, which is the rate of yaw 
sensing device, will pick up components of roll as well as yaw, and at only 
one angle of attack will be a pure yaw sensing instrument as is customarily 


assumed. This fact can be used advantageously to effectively gear the rudder 
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to the ailerons by making the yawing moment due to roll positive. This elim- 
inates adverse yaw upon entry into turns and rolling maneuvers. A stabilizing 
device can be desimed then that will improve the damping in yaw to meet the 
lateral oscillation requirements, making the aircraft spirally stable at the 
same time, and also improve the response to the rolling control. The problem 
of maintaining equal effectiveness of this control with variations in speed 


is also discussed in this analysis. 
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DESCRIPTION OF THE AIRPLAN® AND THE DAMPER SYSTEM 


The airplane under consideration is a single seat, midwing, high 
performance jet propelled fighter with a conventional control system It 
has a two piece rudder actuated by cables from the rudder pedals, not power 
boosted. The physical characteristics, dimensions, and stability derivatives 
were supplied by the manufacturer and checked with independent wind tunnel 
data or methods of calculation as specified by the NACA. 

The damper system to be applied in this case would be known as a 
force type as opposed to a power-boosted type. This means that the servo 
drive has to have an output of the same order of magnitude as a human pilot's 
control effort. In the other case, the servo drive can operate into an arti- 
ficial feel system, boost, or irreversible control power control system Of 
course, the time lag of this force type system is not increased by the in- 
clusion of the time lag inherent in hydraulic rudder boost systems. Consid- 
erable range in power size is available to meet varying requirements. Time 
lags for a force type system can be reduced to less than one tenth of a second. 

The proposed damper system is comprised of three principal parts. 
(1) A rate gyro, a single degree of freedom gyro mounted with its spin axis 
in the XZ plane of the aircraft at some angle to the body reference axis. 
It is constrained by a spring control whose restraint is caused to vary to 
compensate for speed changes. The gyro produces a precession displacement 
and a corresponding electrical signal proportional to the angular rotation 
of the saat in space. (2) A controller-amplifier section, which receives 
the rate signals from the rate gyro, amplifies this signal and imparts phase 


f 


lead. It contains a derivative network which will not pass steady signals 
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from the rate gyro. This network prevents the damper froin opposing the 
pilot in a steady state maneuver such as a turn or roll. Any transient dis- 
turbances in the steady state maneuver will be damped out, however, and the 
beneficial damper contribution upon entry into these maneuvers passes. The 
gain and the phase lead are adjustable to obtain the optimum settings for 
particular applications. The compensating spring restraint feature could 

be included in this unit instead of at the gyro. A typical unit would also 
house the rate gyro in its chassis. (3) The servo drive supplies the power 
to move the rudder upon signal from the amplifier-controller. If used in an 
irreversible application (which it is not in this case) the feedback from 
the load, or position follow-up would not be required since the artificial 
feel spring systems will recenter the load when the signal vanishes. 

The total weight of this system could be kept to within fifteen 
pounds, and the space required to less than one third of a cubic foot. The 
gyro should be located near the center of gravity so as to escape the effects 
of body bending. In this case body bending should be a negligible. This 
problem is treated in reference 8. 

It is appropriate that the distinction between an ae and a 
stabilizing device should be mentioned. The autopilot has authority over 
the aircraft, presumably by the choice of the pilot. He may overpower it or 
control it through autopilot control settings, but the autopilot has direct 
control of the aircraft. A stabilizing device, such as the proposed yaw 
damper, operates in conjunction with the pilots the pilot, however, has con- 
trol and is superimposing his control motions over the damper. The damper 


can be thought of in the same sense as the vertical tail in that it is a part 





PRINCETON UNIVERSITY PAGE 
AERONAUTICAL ENGINEERING LABORATORY . REPORT 172 


eS tp A RE Lp : 





a te cer eee ee ee St 
re 








of the aircraft which is in continual operation contributing stability. It 
is the opinion of the author that it should be an absolute requirement that 

a stabilizer shall have no feedback to the control system of the pvilot; the 
pilot should not be conscious of the damper's operation except that he feels 
that. the flying qualities are improved, and that he is able to fly a more 
highly damped and smoother path. This isolation can be easily accomplished 
in an irreversible control system. In the application under consideration 

it poses a difficult problem that can probably be best solved by dividing the 
two piece rudder between the pilot and the damper. That is, using a control 
surface for the damper alone. No solution in which the pilot can feel the 


damper in operation through the rudder pedals should be considered satisfactory. 
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LATERAL DYNAMIC CHARACTERISTICS 


The first step in the problem is to determine what the ealeulaten 
lateral-directional oscillation characteristics of the aircraft are. 

Three configurations of the aircraft have been considered. 

I. The maximum cruise configuration, at 35,000 feet, clean, at a Mach 

number of .68. 

II. The 50% Normal Rated Power configuration at 10,000 feet, clean, at 

a Mach number of .3c. 

III. The Power Approach condition with flaps and gear down and with power 
for level flight at a reasonable approach speed. 

These three conditions give a sufficient range of altitude and speed 
to reasonably judge the lateral dynamics of the aircraft with and without a 
damper system. Configuration II was selected because pilots reports have 
indicated that at that particular altitude, at Mach numbers below .6, the 
Dutch Roll or snaking has been particularly noticeable. Conditions I and 
Iii represent the extreme opposite ends of the speed and altitude range in 
configurations of operational interest. 

The equations of motion as presented in reference 9 have been used, 
with the exception that the effects of the product of inertia have been in- 
cluded only in condition I. In this condition ™% is negative... Im the other 
two conditions | is positive and it is conservative to neglect the product 
of inertia. For this particular aircraft the product of inertia was found 
to have little effect. These calculations are reviewed in the section on 
calculations. The data is presented in Table II and plotted in figure <. 


All the points lie in the unsatisfactory region. Flight tests have 
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shown good correlation with these points. The calculations Bae’ comlirol 
fixed, iniicating an unsatisfactory Dutch Roll characteristic. The boundary 
between the satisfactory and unsatisfactory regions of figure 1 of reference 
l is currently under at iigattion. The calculated voints show the aircraft 
to be in a region which has been called tolerable, meaning tiring or un- 
pleasant, but not necessarily dangerous. The calculations show that a con- 
dition exists which would certainly be annoying and detrimental to the 
earrying out of the airplane's mission. 

Whether this condition is of sufficient severity to warrant a change 
or modification in already operational units is doubtful. That is a question 
of opinion of the pilots operating the airplane. All indications are that 
future designs are going to need automatic stabilization; analyses of this 
type will be called for. Applications to longitudinal motions are also in 
use. Future airplanes may well use automatic stabilization both laterally 


and longitudinally with the pilot imposing his control efforts over the 


automatic controls. 





=. 


PRINCETON UNIVERSITY PAGE /O 
AERONAUTICAL ENGINEERING LABORATORY REPORT (72 


NR I: RD Pa A 


CONTRIBUTION OF THE YAW DAMPER 


Referrinz to figure lb the relation between the body reference 
axis, the gyro spin axis, and the wind or stability axis is shown. Sub- 
scripts R (reference), G (gyro), and V (stability) are used in the following 
development. 

For simplification the time lag of the yaw damper system is first 
assumed to be zero. It will be later shown that the critical time lag of 
the aircraft-stabilizer system is considerably higher than the time lag of 
any mechanism under consideration. 

The equations of motion are written with respect to the wind or 
Stability axes. The derivatives contributed by the damper must therefore 
be related to the V axis for inclusion in the equations. The rudder, how- 
ever, introduces yawing and rolling moments about the reference or body axis 
R. Also the yawing and rolling moments so introduced are proportional to 
the rate of yaw as measured with respect to the gyro axis G. 


Then the development of the contributions of the damper is as follows: 


ad 
a 
K = ’ rad/\ ad /sec control gearing or gain. 
(Cp Sue) Set Ce hy 
but 


Y * Yay cos® + A, sin § 


and 
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C. = (Cy = (Cs cos a+ (Cy ys sin a 
Cy = (Co )y = (Cy )e cos a + (Cd sin a 
substituting into the above equations, 


BASF el Foye 9 +f, sing ] 

+K eS ane +f, sin @ ] 

Gy =K Cp [vy cos f +f, sin & ] 
r 

My [Y ,cos$ +f, sin¥ ] 


z 


These expressions can be seen to have the form of to 


en : dc ; 
4c = iv AY + ayo? 
aC dC, , 


ear |e + ag 42 
therefore 
SC =K cos iG + 6 a} 
ny E n ; "5 
-_ =Ksin@ [C+ By. a] 
"r r 


AC), =Kcos€ [C, + Ca] 
ky < L< n 


xz r 
a, =K sing Gy .* Ri a] 
Tr cr 


since 


PAGE i" 
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cos a 


cos @ 


sin a 


tal derivatives, where 
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ny ny, =! n.. LV 
2V 
- Av 
e > K cos & “a 
r 
a 
Ac 5 K gin § “ns 
P r 


neglecting the product aC » Which is small compared to C : 
R Nk 


r r 
The damper contributions are seen to be directly provortional to 


true airspeed. The rudder power also is seen to vary with altitude from 
Table II. 

The rudder power coefficient has a negative sign, and the rolling 
coefficient due to rudder deflection has a positive signe It can be seen, 


therefore that if 4 has a negative sign the increments Me and acne will 


Y 
desirable sims for these contributions by the damper. It is therefore 


be positive, and the increments ee and M“,. negative. These are the 
concluded that the gyro axis should be so set that it is always below the 
stability axis, ( 8 a negative angle) as shown in figure lc. 

The increnents Mo. and ami have little effect on the lateral 


dynamic calculations, and sum Cy + ny a} is considerably less than 


r 
[c., + Cy a]; for these reasons the conservative assumption is made in 


r r 
further calculations that only AC. and “yy supplied by the yaw 
a 


damper. 


The increments added by the yaw damper are then assumed to be 
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A «= & cos 3 C 
nz n 


Y Oo. 


a) 


AC. = Ksin€ C 
ny x n 

r 

The optimum settings of K and 8 =f(a,D) are discussed in a 


later section. 
Reference 10 treats the effect of C_ (= Oa on lateral dynamic 
p 


stability. Making C_. positive was shown to have a beneficial effect in 


ne 
shifting the points on a plot such as figure 4 to the satisfactory area. 
As K and -& are increased, the positive increment of Ce. lengthens the 
period at substantially the same time to damp to one-half amplitude, improv- 
ing the cycles to damp to one-half amplitude. This increment in oo nust 
be limited, however, since too great an amount leads to a long eid un- 
stable oscillation of the roots that were identified as the rolling converg- 
ence and the spiral mode. These roots first become a long period stable 
oscillation and then as -. is further increased the oscillation becomes 
unstable. This shows that the product VK sin § must be limited so as to be 
satisfactory in configurations as contrasting as I and III. 

The reasons for seeking me contributions from the damper are to 
assist in correcting the Dutch Roll oscillations and to oppose the adverse 
yaw introduced by the ailerons and the normal negative derivative (adverse) 
of yaw due to roll. By overcoming the adverse yaw in this manner, entry 
into rolls and turning flight can be considerably improved. This benefit 
can only be enjoyed on the entry into these maneuvers, since the steady yawing 


velocity signal will be intentionally blocked to keep the yaw damper from 


opposing a steady turn initiated by the pilot. 


aa 
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ee is seen to have its greatest influence in the by term of the 
stability eae As this derivative becomes more positive it drives by 
to zero. As this occurs the rolling and spiral modes first become a stable 
oscillation, and finally a long period wstable oscillation. The limiting 


value of -§ for the fixed K was then calculated by setting b, to 


1 
zero, keeping all derivatives except Ow, and Dan, the same, and introduc- 
ing these as £(§ ). The limiting P determined in condition I is shown 
in figure 5. Since a varies over a small known range, and g z=a-LT, 
it is possible to estimate the range of effective g - See figures lc, 3. 
AC. and x » the non-dimensional forms of the damper contribu- 
tions, menveaite to vary with speed. Therefore in order to maintain equal 
effectiveness of the stabilizing control throughout the speed range, K 
must be made to vary inversely with speed. One convenient way to do this 


Le to the spring restraint of the gyro, 


is to apply the inverse of gq 
thereby modifying the electrical signal from it. The V's will cancel, but 

a varying factor of ole remains that will have to be calibrated out as 
well as possible. A close approximation to the desirable variations in K 
with the three conditions is shown in Table III (K to make C1 fo" 05). Com- 
paring these values with those shown in Table II (calculated by using a 

gn i/2 variation) it is seen that the contributions tend to be too large at 
altitude. Figures 2, 5, and 4 confirm this. Applying dynamic pressure from 
a pitot source seems to be the most convenient way of compensating, however. 
A non-linear spring and bellows with a pitot tube source such as are used 

in artificial feel systems would probably be the best solution. This penalty 
in gadgetry must be paid to realize the benefits of the positive yaw due to 
roll introduced by the damper. Positive Cs must be limited, and the goal 


p 
should be similar performance at all altitudes and speeds. 
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SETTING GYRO AXIS ANGLE AND GAIN 


The question of just how damped,a motion the pilot would like to 
have is still under considerable discussion. There is no question that in 
a system of this sort it would be possible to add enough damping with a 
Stabilizer to the aerodynamic damping to critically damp all lateral oscil- 


lations. However, the standards for the performance and stability of the 


4 
aircraft-stabilizer servomechanism are not as yet established. See refer- 


ences 11 and le. 

This analysis was made on the assumption that the gain should be 
the minimum required to bring the cycles to one-half amplitude to one half. 
There is evidence to show that the motion will then appear critically damped 
to the pilot. Reference 8. 

A Routh's discriminate investigation using K and as variables 
is possible, but exceedingly tedious. Determination of the optimum values 
of these variables by a plot of Cif vs VD with K as a parameter is 
shown in figure 3. This approach is relatively easy and equally informative. 

Referring to figure 3, curves for K =1 and K =@ _ in condition 
Iii are shown. The higher value of K does not give sufficiently better 
results to warrant the higher gain settings. It is in this condition (sea 


level at approximately 1.20 V ) that the basic value of the gain, or 


stall 
‘control gearing, is established. At higher speeds (or q's) the gain is 
reduced as already explained. 

Selecting the basic value of K = 1 in condition III, the effective 


value in condition I becomes .5, in condition II, .817. The curve for con- 


dition I for this effective K is shown up to the limiting value of [PD , 
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determined as discussed in the previous section. In this particular case 
there is no limit to the gyro angle in condition III. In other words with 
gains below 2, too much +C cannot be introduced. The effect of holding 
K constant and inclining the gyro axis to the point that damper conti imax 
tion is entirely aC. is shown in figure 4. The points of g =O sand 
g = -90° on each curve contrast the effects of stabilizers sensitive to 
pure yaw or pure roll. 

Selecting a gyro angle of 30° and K = 1, the range of values of 
Ci for a variation of angle of attack of fifteen degrees is shown. This 
setting of D and K is taken for further investigation since it is the 
maximum practical gyro angle at this K that will not introduce an unstable 
oscillation due to too much +0 in condition I. Some lower angle might 
prove more desirable from Plight tests. This is a question of how much 
positive yaw due to roll is best for entry into rolls or turns. From figures 
8, 9, and 10 it can be seen that the response to the aileron control is 
radically modified. 

In general, a quick procedure for selecting gyro angle and gain is 
as follows: | 

(1) Using the single degree of freedom equation in yaw, a close 
estimate of the total amount of damping required to damp to any given vi /2 
is easily obtainable. From this, the amount of sony that is to be added 
by the damper can be estimated. Knowing the rudder power, the approximate 
gain to get this increment of damping can be found. 

(2) From the adverse yaw and yaw due to roll coefficients the total 


amount of ~~ can be estimated. Using the gain calculated above. 
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the required E to get the increment in me estimated above can be solved 
for. using the expression obtained under the damper contribution section. 
For the gyro angle, LD =a - ¢ . 

The compensating device on the gyro will keep these increments essen- 
tially constant with speed and altitude. To check to see if an instability 
has been introduced, select a high speed condition where the final ratio of 

ng °° cons is high (usually at high altitude and speed) and calculate the 
by term of the stability quartic. 

Final checks using the three degree of freedom equations and involv- 
ing a separate solution of the quartic for each point as was done in i cme 
5 can then be made for a critical configuration. The assumption of zero 
time lag will have to be kept in mind and later checked. At this point it 
probably would be cesirable to begin flight testing the system for the final 
selection of these paramenters. 

In order to prevent any feedback through the pilots rudder control 
system, the yaw damper for this application will probably have to have its 
own rudder area. The calculations in this report have been made on the 
assumption that the whole rudder was available for the yaw damper. This 


can be done when the movements of the rudder caused by the damper can be 


masked from the pilot through a boosted or artificial feel control system. 


Spf, 
ewer + 
ns Ls Sb 
YP a. 
and 
Seh 
C D -” h 


= C ——=.—(C 
ks a? I, 2S 
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The terms involving Cy have been neglected in the development 


- 
of the contribution of the yaw damper. Since 


AC = K cos f C 
ni n 
- 6. 


A = Ksin€@ C 


C e 
ms ne 


rE 


the same results can be obtained by keeping the product of KC constant: 


r 
or as the area Sp° for the damper, is reduced, the gain K is increased. 


Note that the vertical location of the area for the damper is a negligible 
effect. The maximum practical value of K obtainable with this type of in- 
stallation (probably 2 or3)will limit the nininun size of the damper rudder 
area. 

In summary, the optimum gyro angle is chosen as 50 depraie, and the 
optimum gain as 1 at sea level and approach speed. This settine has the 
following results: 

(1) Moves the points shown in figure 2 from the unsatisfactory 


region to the satisfactory damping region where C is less than .5. 


We 

(2) Adds a sufficient amount of = NG in each case to make the 
aircraft spirally stable in all conditions. 

(3) Adds a +C, increment of an approximately constant value to 
ail conditions, which is considerably larger than the total adverse yaw due 
to roll and the rolling control. 

The beneficial effect of the additional damping is shown in the 
response data plotted in figures 5, 6, and 7. The improved response to the 


rolling control is shown in the response data plotted in figures 8, 9, and 
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10. 

The effect of varying VY at a constant K for condition III is 
shown in figure 4. 

Figure 3 shows the selected points in the three conditions, and a 
Cie vs bP plot for constant K for conditions I and Lil including the 
range of angles of attack. 

All calculations involving the contributions from the yaw damper 
were made without using the small angle assumptions or neglecting the influ- 
ence of roll due to rudder deflection. The simplifying assumptions made in 


the analysis of the damper contributions were found to be justified. 
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It was pointed out in the Introduction that aerodynamic damping 


should be augmented by gome other means than increasing the tail length or 


area. Increased tail length is limited by operational factors such as car- 


rier elevator sizes for naval aircraft. The amount of vertical tail area 


also has a practical limitation which precludes the increases in damping 


that are desired. 


The yaw camper has the advantage of increasing C, without affect- 


Yr 


ing the value of Co - Tail size can then be chosen without consideration 


of the requirement of high damping, and . can be set from consideration 


of static stability, stalling qualities, etc. Dismissing the possibility 


of adding to conventional tail length, the vertical tail area would have to 


be increased by a factor of about four to achieve the same results in damp- 


ing that the addition of this stabilizer will. 


Reference 8 points out that with the XB-47 it was physically impos- 


Sible to change the effective dihedral negatively enough to correct the 


lateral dynamics because of ground clearance, and that the tail area would 


have to be at least doubled to add the necessary damping. It was 


explained 


that in general, while increasing tail size increases damping, the gain in 


smooth flying qualities is questionable. While the tail size increase adds 


\ 
to the damping, the excitation of the aircraft due to side gusts increases 


in the same proportion. 


This XB-47 report uses the servomechanisms approach to stability 


analysis as outlined in references 11 and 1. In one polar stability plot 


White shows that the effect of increasing C and C. together (as in the 


Yr 
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case of increasing tail area) results in little or no improvement in damping 
because the resonant frequency of the aircraft has been increased. 
Consider the single degree of freedom equation in yaw, 
a = a? ee ee Oe 
B 0 
r 
This equation may be transformed into the standard servomechanisms 


form using the following definitions: 








seh 
damping ratio 4 4 a 
20 
> 
undamped natural frequency uf = +t 
2 
ns 
rudder effectiveness i C ‘ 
"B 


The equation of motion then hecomes 


Yeh oY +¢u% = Ku S 
O O s “0 
and the 


damped natural frequency at L1-% ; 


resonant frequency wo, = @ f1- 298 


For the usual values of 4 ‘ Wes Ws and @, may be considered 
equal. White, reference 8, considering the single degree of freedom case, 
sets up his criterion of damoing as the ratio of dynamic yaw response occur- 


ring at w, » to the effective static yaw response considering w= ws 
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The measure of the damping then becomes inversely proportional to the square 
of the undamoed natural frequency, hence inversely proportional to or ’ 

2 “Sa F 
since WA Ww C /i, - 


Oo n 
P 
Also from the single degree of freedom equation. 


5 es V1 - 6? 
1/2 i 

White states that within the limits of engineering accuracy, the 
single degree of freedom equations may be used for calculations involving 
the design of the yaw damper and to obtain Cie ° 


Results from single degree of freedom calculations are included in 


Table III. 
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TRANSIENT RESPONSE INVESTIGATIONS 


In order to graphically present the contrast between the transient 
response of the aircraft with and without the damper system operating, the 
following calculations were made. For condition II the responses +/8. : 
B/E. » and B/E, were found for a disturbance in sideslip effecting each 
of the three equations of motion. The function introduced is a mit step 
function in sideslip applied as a step function at t = 0, and neutralized 
as a step function at t = 1. Condition II was selected because in this 
configuration pilots have frequently reported unpleasant lateral oscilla- 
tions. The forcing function used is an attempt to simulate a disturbance 
that might be initiated by turbulence or by a quick displacement in sideslip 
however introduced. 

The airplane in the original condition is seen to have persistent 
oscillations in all three plots (figures 5, 6, and 7). The period is long 
enough so that the average pilot could damp them out with the proper appli- 
cation of controls; however, requiring the pilot to continually oppose these 
motions is tiring and detrimental to carrying out his mission. There is a 
considerable improvement in yawing velocity and sideslip, but the rolling 
characteristic with the damper in operation should be investigated in flight 
test. Too high a ratio of roll to sideslip is particularly wpleasant to 
the pilot. From data available in other applications, oscillations that damp 
to half amplitude in less than half a cycle appear dead beat to the pilot. 

From this data and the rudder hinge moment parameters, the rudder 
pedal force that must be supplied by the servo power section of the damper 


can be estimated. An estimate of the amount of rudder throw required of the 


PRINCETON UNIVERSITY PAGE 24 
AERONAUTICAL ENGINEERING LABORATORY REPORT /72 


a cee ree cag NN AE SA Po I SO a a eS SG GREENS GS 





damper rudder surface can also be made. The final setting of the stops 
limiting the throw of the damper rudder surface must be determined experi- 
mentally since the control required will depend upon the characteristics of 
the servo system as well as the airplane responses and the frequency range 
to be damoed. 

The response to an aileron step function is presented in figures 
8, 9, and 10. Condition III is investigated since it is probably the most 
unfavorable. Adverse yaw has been included in the calculations. The large 
positive yaw due to roll introduced by the yaw damper is seen to drastically 
change the initial response upon entry into the rolling maneuver. The air- 
plane no longer starts to yaw in the wrong direction for the turn, and the 
roll builds up rapidly. The sideslip builds up at a slightly faster rate 
also, but not rapidly enough to offset the good effects in yaw and roll. 
The derivative network in the yaw damper circuit will have an RC circuit 
with a time constant so desiged as to only pass the first portion of a second 
of the contribution of the damper so that it will have the effect of assist- 
ing in the rapid build-up of roll without the adverse initial yaw, and still 
not oppose the steady-state maneuver. In this manner the damper should con- 
siderably improve the response to the rolling control. 

From figure 5, (the response in yawing velocity to the disturbance 
in sideslip) the maximum yawing velocity per degree disturbance is seen to 
be approximately two degrees per second per degree sideslip. Then ten degrees 


of rudder would be required for a gain of one and an initial disturbance of 


five degrees, since 6 e lao /a | Y 
D 


Considering a ten degree throw in either direction as the limit for 
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damper applied rudder deflection, it can be seen that the "with damper" 
curves shown in the aileron response plots will be cut off by the rudder 
reaching its stops. This cut-off point will be determined by the amount of 
aileron applied since figure 9 is presented per degree deflection. When 

the product of the gain and the yawing velocity exceeds the available full 
deflection, the rudder will be initially held against the stop with a ten 
degree favorable deflection. The RC circuit will then cut out the damper 
rudder contribution as the steady signal is sensed by the gyro. Should 

there be a tendency to oscillate, the servomechanism will damp the variations 
about the steady maneuver out. 

The curves showing the improved performance with the damper oper- 
ating (figures 8, 9, and 10) will then be modified by gervo rudder stop 
setting. Startling improvement in the entry into rolling maneuvers should 
still be obtainable with limited rudder throw. 

Two factors should be here reiterated. First, these calculations 
were made for the maximum value of gyro axis inclination, and hence positive 
C. increment, possible at the selected value of gain. In other words these 
een responses are for the maximum gearing of the rudder to the ailerons. 
Some smaller angle of gyro axis inclination might prove more desirable from 
flight test or more exhaustive calculations. Second, the device introduced 
to assure nearly constant oe ae the damper with variations in 
speed and altitude will automatically reduce the gain at favor sciseities. 


This should prevent destructive rudder motions as it hits its stop after the 


gyro senses a large rolling velocity. 
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DETERMINATION OF CRITICAL TIME LAG 


In all the calculations to this point the assumption has been 
made that the time lag of the damper system is zero. The critical time lag 
of the aircraft-autopilot combination is now investigated to determine how 
much of a limitation this assumption imposes. 

The frequency response method outlined in reference 14 is applied. 
The critical time lag of the system is defined as the time lag that results 
in a neutrally stable or steady~state oscillation. The steady~state response 
Ky = IY /5,| of the airplane to a sinusoidal forcing function of unit ampli- 
tude was first calculated. The calculations were made for condition III 
since in this condition the ratio of the damping to the critical damping 
was lowest. The gain of the yaw damper is known, having been set at K=l 
for this condition. The gain, K. is the amplitude of the control-surface 
oscillation produced by the yaw damper mechanism in response to the oscil- 
lation of the airplane yawing velocity. Intersecting K | = [1K | with the 
airplane gain or amplitude plot the maximum permissible time lag in the 
system can be determined. The slone of the line through the point where the 
phase angles are equal is a measure of this critical time lag. The condi- 
tion for stability is that at each angular frequency where the phase angles 
are equal, the ratio of the amplitudes of the airplane to the autopilot must 
be less than one. 

From figure 11 the critical time lag of the system chosen is seen 
to be 1.ek seconds. Critical time lags for gains of two and three are shown 
also. For a damper gain of three, the limiting time lag is .7 considerably 


in excess of that inherent in typical force type systems. The critical lag 


&é 
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an refex ance 8 is a little over .5 seconis for K =1. For servo- 
chanisms with larger time lags introduced by boost systems, this calcu- 
on would undoubtedly show a closer margin. In this case, the time lag 


tem is not a Limiting consideration. 
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As a result of the theoretical investigation of the application 
of a yaw damping servomechanism to the aircraft under consideration, the 
MT ing conelusions are reached: 

(1) This plane does not satisfactorily meet the requirements of 
reference 1 in regard to lateral-directional oscillations. Control-fixed 
calculations predict that the airplane will exhibit a "Dutch Roll" in all 
the conditions investigated to a degree that could be classified as objec- 
tionable but not dangerous. These calculations substantiate pilots reports 
of “snaking” or "Dutch Roll." 

(2) By installing the yaw damper system as proposed, the aircraft 
will then meet the requirements of reference 1. The lateral oscillations 
Will appear nearly critically damped to the pilot, and he will be able to 
carry out the mission of the aircraft more effectively with less effort. 

(3) The damper will siepealt enough damping in yaw to make the 
aircraft always spirally stable. It will also introduce a constant positive 
yawing moment due to roll which will improve the turning and rolling per- 
formance. The system will be so designed as to improve entry into steady 
turning or rolling maneuvers, but not oppose them. 

(4) Damping introduced by a servomechanism which applies rudder 
control proportional to yawing velocity is probably the most promising if 
not the only satisfactory way to augment the aerodynamic damping. This sta- 
bilizer introduces damping without increasing the resonant frequency (or 
the effective static yaw stability), which insures results not obtainable 


by increasing vertical tail area. 


ccc 
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(5) For this particular installation the optimum settings of gyro 
spin axis angle and damper gain are Db = 30° and K=1. These values have 
been checked in a limiting condition and the amount of positive yaw due to 
roll introduced is within the limit that would cause a lonz period unstable 
oscillation. Some lower value of UD’ that would still produce a +O might 
prove more desirable from flight tests of the effect on rolling maneuvers. 
In all cases, the spin axis of the gyro should be kept below the stability 
or wind axis. Figure lc. 

(6) It should be required in the case of a stabilizing system 
(as opposed to an autopilot) that no feedback to the pilots controls be 
allowed. In this case the most practical solution appears to be the use of 
a separate rudder area for the yaw damper mechanism, probably splitting the 
two piece rudder between the pilot and the damper. A differential linkage 
or gearing, with the pilot and stabilizer as independent inputs with pro- 
vision for balancing out the servo applied hinge moments could be developed. 

(7) Calculations show that the damper system as proposed will 
improve the transient response to a disturbance in sideslip, and also improve 
the response to an aileron step function input. 

(8) The critical time lag of the stabilizer-airplane system is 
considerably in excess of the maximum time lag to bea expected of any servo- 
mechanism of the type vroposed. The assumption of zero time lag in the cal- 
culations proved to be reasonable. 

White, in reference 8, has an excellent presentation for the invest- 
igation of the stability of the stabilizer-airplane combination which he calls 


the inverse frequency response method. It is felt that this type of approach 
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should be used after the experimental frequency response data of the actual 
servomechanism to be used is available. Otherwise results would be of 
academic interest only. The method presented here will enable the designer 
to quickly select a system with:the required gain, time lag, and power out- 
put from the servo drive and install it near the optimum position relative 
to the airplane reference axis. Flight testing the system from that point 


would probably yield the quickest dependable results unless an analog com- 


puter were available. 
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- CALCULATIONS 


The lateral equations of motion are given in reference 9. In all 
cases Cy : Cy - and XY were assumed td be zero. The equations including 
r 
the product of inertia were used only in case I. Otherwise they were used 


as follows: 


it 
O 


(m_D - y,)p - C2 + md? 


fi 
© 


“Cy ,B (1,D*- oP) ¢ GD 


-C, fe ODP GL Uke DY weve 


B 1 
then, 
b ps web » See b N+b. 2A=O 
4 by 2 1 oO 
where 
by - att, 


b “(6% 5 + “yg x! : Mon 


bo =C, Cy 1 + (mp, + Cy LC +mic -mo,.c, 

2 “Ye oe O b ee ny os iy 0 4y ng 

oe ei-C. C, C +mcC,C = pwee.© @~- C_C, I + C.C C 
Ion, Le Q Y 4G) Tm 

i Ye g e*. oO L2& “2 p y 2g 

b @C.(C, C_. ~Cy. C ) 

O L Ke ny, fy Re 


The fastest and most accurate method for the solution of the stability 
quartic was to use a combination of the Lin approximation method and the 
Zimmerman plot as explained in references 11 and 13. With the Lin methad, 
values could be obtained that would lead to an intersection in two or three 


points with the Zimmerman plot. Im general, there were two real roots and a 


eee 
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complex pair, although where a large value of TC was used there were 
two complex pairs. 
Where the roots were Ny h and a+ if, 


_ 2693 
tie a. 





seconds 
Pp = ats seconds 
B 
. p 


For the plots in figures 3 and 4, where each point represents the 

solution of a quartic, the equations for b, 3.2.1.9 were used where the 
29 y= 'S 

variables were C.- and ony determined from the particular value of 
5 = é fa.) for™Gaat point. 

Using the above three equations of motion, a general solution for 
forcing function inputs of A. B, and C respectively was made. The transient 
responses took the following forms: 


A fs) 2 
Wy ag” + an A” + a, /\ + ie 








6(t) A 
where 
"5° p's 
- ag) J obec. - C(mc, + Cyl) 
Rg : : XY 
ne x O ng O he x x 
a, = A(Cy Cow = oe Co. )-BC_ Cy + C09, Cy 
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In the solution for the disturbance in f, dt) became B , and A, 


B, and C were CL, CO and C. 
¥ co” n 
B p B 


In the response to the aileron step function, SH) became a and 


A, B, and C were Oe Co 3 and C ® 
6 Oe 


& & 


4, and hence the roots ” Ns » and a F if are Inown. 


Ll 9 
The solution for the typical lateral transient response takes the 


form 


Peat 
Ect ye © aC Ff c, ¢\? + Bp° sin( Bt + F) 


12 
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By use of the Heaviside Expansion Theoren, 


a 
C = 7 
oO 

O- 


Z, e 
6 te as i. i a, i. a, ~~, 
r 3 2 
Dr (4b, - 3b, AY + 2b, ie +a) 


Then from end conditions 





- |e + C, + C,) 
gin F 
es 5: -B(C. + Cy + C,) 


C. =o, ~ Nace HC, + Ge 


where C. represents the initial condition of Y op, OF g at t #6. 


C. = +~ . An initially relaxed system is assumed. 


To find the frequency response for the critical time lag calculations, 
standard servomechanisms methods were used. It is believed that the best 
method of reducing the data is outlined in reference 8, and is as follows: 

For the airplane, 

. 3 2 
ye 2 aD + aod + a,D + ay 
o 4 3 2 
x b,D + b, D + byD + b,D + bo 
Using the general equations developed previously with forcing func- 


tions A, B, and C as Cy gin wt, Co sin wt , and C sin ot the a's 
S é n s 
r 


r r 
can be quickly calculated. 


- Then Ky » the steady-state frequency response of the airplane, can 
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be expressed 
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where the substitution D = iw is made. 


a 
~ Bo + a 


x 


1 oO 


> 
XY) “a, a ao 


: 19,-P2)0 


" 


, 


us ~ bw" + b 
O 


) 


5 
~b,W + by 


2 
Xo 


ain * 


Then for each frequency the magnitude of the airplane gain or ampli- 








tude is 
W x, cos B. 
LA = oo Xo cos B, 


and the phase angle is 


B= B, - Bo 


10. 


ll. 


12. 


15. 


14. 
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PRINCETON UNIVERSITY 


AERONAUTICAL ENGINEERING LABORATORY 


Principal Dimensions 


S = 250 sq. ft. 
b = S5.e0 ft. 

¢ = 7.45 ft. 

pl = 18.45 ft. 

S = 36 sqft. 

S = 8.54 sq. ft. 
Cc = 7o Tt. 


b e 742 SQe fhe 





WI = 12,009 pounds plus fuel 


Total control throw for the rudder, 25°4 1° each Way 
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